Diffusive behavior of a greedy traveling salesman 
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Using Monte Carlo simulations we examine the diffusive properties of the greedy algorithm in 
the ci-dimensional traveling salesman problem. Our results show that for d — 3 and 4 the average 
squared distance from the origin (r 2 ) is proportional to the number of steps t. In the d = 2 case 
such a scaling is modified with some logarithmic corrections, which might suggest that d — 2 is 
the critical dimension of the problem. The distribution of lengths also shows marked differences 
between d — 2 and d > 2 versions. A simple strategy adopted by the salesman might resemble 
strategies chosen by some foraging and hunting animals, for which anomalous diffusive behavior has 
recently been reported and interpreted in terms of Levy flights. Our results suggest that broad and 
Levy-like distributions in such systems might appear due to dimension- dependent properties of a 
search space. 
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I. INTRODUCTION 

Although the classical random walk problem is rather 
well understood by now [l], some of its modifications, re- 
lated, for example, to disordered systems, are less clear 
. Departure from the classical random walk manifests 
often as an anomalous diffusion that appears whenever 
the typical square displacement of a tracer particle R 2 (t) 
does not obey the normal diffusion law R 2 (t) ~ t a with 
a = 1. Deviations include numerous examples of subdif- 
fusion (a < 1) such aSjfor example, photoconductivity of 
amorphous materials Q, diffusion in convective rolls 
or a random walk on the percolation cluster On the 
other hand, there is much rarer superdiffusive (a > 1) 
behavior, of which the Richardson diffusion in turbulent 
fluids is a classical example @ . Some macromolecules or 
biological systems also exhibit anomalous dynamics with 
superdiffusive behavior, as reported, e.g., for micelles Q 
or living cell migrations [8[. 

Recently, a number of studies have indicated that for- 
aging or hunting strategies of some animals also show su- 
perdiffusive behavior and examples include ants Q , fruit 
flies [loj ]. or wandering albatrosses [ll|. The prevailing 
understanding is that while searching for food these ani- 
mals move in steps of varying length with fat power-law 
distribution consistent with the Levy flights. However, 
some other works show that this hypothesis should be 
taken with care. Indeed, re-analyzing the data for the 
wandering albatrosses, Edwards et al. concluded that the 
previously suggested Levy-flight scenario appeared due 
to a spurious inclusion of very long flights [12]. It has 
also been claimed that the idea of animals using stochas- 
tic but essentially blind search strategies does not seem to 
be viable since it neglects the role of animals' intelligence 
and experience in guiding them [13| . Nevertheless, some 
more extensive recent studies for 14 species of open-ocean 
fish of prey seem to support the Levy-flight scenario [3] . 

In the present article we show that a simple search 
strategy, i.e.: go to the nearest (intelligence) and not yet 
visited (experience) place, leads to an interesting diffu- 



sive behavior. In particular, a normal diffusion that takes 
place for d — 3 and 4 systems is modified by some loga- 
rithmic corrections in the d = 2 case, which suggests that 
this is the critical dimension. The distribution of dis- 
tance does not have fat power-law tails, however, when 
the search exploits all available places, much broader and 
possibly Levy-like distributions can be found. 



II. THE MODEL AND THE NUMERICAL 
METHOD 

We consider a traveling salesman [ll| that has to 
visit cities which are uniformly distributed in the 
d-dimensional space of linear size L. We assume a unit 
density of cities, i.e., N/L d = 1, where N is the num- 
ber of cities. The salesman is not interested in finding 
the minimal tour but visits cities according to the fol- 
lowing greedy algorithm: in each step, from the last vis- 
ited city it moves to the nearest city which has not been 
visited yet. For most simulations, the starting point is 
the city placed in the center of the space, i.e., of the 
Cartesian coordinates (L/2, L/2, L/2). The quantity 
of main interest is the average square distance (r 2 ) from 
the starting point and especially its asymptotic form af- 
ter many steps t. An example of a typical trajectory 
in the two-dimensional case is shown in Fig. [T] Let us 
emphasize that for a given starting point and a distribu- 
tion of cities, the route of the salesman in our model is 
completely deterministic even though it resembles a ran- 
dom walk trajectory. Our model might be considered a 
version of the so-called deterministic random walk [161 ] . 

To examine the asymptotic behavior of (r 2 ), it is es- 
sential to average over many configurations with large 
number of cities N so that a relatively large part of the 
trajectory of the salesman was not affected by the bound- 
aries. Because of these finite-size effects, we terminate 
the generation of trajectories after number of steps that 
is much smaller than N. As a result, the t-th city on a 
trajectory is typically a close neighbour of the (i-l)-th 




FIG. 1. An exemplary configuration of the path of the greedy 
traveling salesman for N = 2 x 10 4 cities (dots) and duration 
t = 10 3 . 

city. Instead of searching for the nearest not yet vis- 
ited city in the entire configuration, it would be much 
faster to restrict the search to the cities which are in 
the vicinity of the last visited city. Such a modification 
can be implemented by means of the sectorization of the 
space. In addition to coordinates, we store the informa- 
tion about a sector to which a given city belongs. Then, 
to find the t-th city, we search the sectors surrounding 
the (t-l)-th city (see Fig. [2j. Let us emphasize that the 
above described sectorization does not introduce any ap- 
proximation as for a given configuration of cities the al- 
gorithm generates exactly the same trajectory as the one 
without sectorization. In our simulations we used such 
mesh that each sector contained approximately 10 cities. 
(In the two-dimensional case, where the longest trajec- 
tories had to be generated, the sectorization reduced the 
simulation time by the factor ~ 10 3 .) 

III. RESULTS 

Using the above algorithm, we simulated d=2, 3, and 4 
systems containing up to 3 x 10 s cities. Results for the 
two-dimensional case are shown in Fig. [3] Since the 
asymptotic slope seems to be bigger than unity, these 
data show a slightly superdiffusive behavior, but more 
quantitative analysis of the slope would be desirable. 
One should take into account that the asymptotic regime 
in the numerical data will be visible only in an inter- 
mediate interval since both the initial transit as well 
as large-t data should be neglected (for long trajecto- 
ries (r 2 ) is strongly affected by finite-size effects). To 
estimate a, we first examine the time-dependent quan- 
tity a(t) defined as the local exponent estimated from 
data spanning over one decade and centered at t (i.e., 
to estimate a(t) we used (r 2 ) calculated for t' such that 
log 10 (t) - 0.5 < log 10 (i') < log 10 (t) + 0.5). 

The obtained results are presented in Fig. 0J As ex- 



(a) (b) 

FIG. 2. Sectorization of the two-dimensional space. For a 
trajectory that at a certain moment arrives at the sector so, 
the program searches for available cities within so- (i) If it 
finds some, then to find the nearest city it is sufficient to 
examine so and its neighbouring sectors as indicated in (a). 
This is because any city outside these sectors is further from 
any point in so than the length of the diagonal of so (which 
is the largest distance within so), (ii) If there are no available 
cities in so, the program searches for them in a larger area 
(the bold square in (b)). If it finds some, then to find the 
nearest city it is sufficient to examine sectors as indicated in 
(b). (iii) If no cities are found, the program examines even 
larger areas (until some cities are found). 
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FIG. 3. The logarithmic plot of the time dependence of (r ) 
for several values of N indicated in the right bottom. The 
straight line corresponds to the linear increase (r 2 ) ~ t. For 
each N, the plotted results are averages over at least 10 
independent configurations. 



pected, for large t all curves bent downward and in the 
middle range a plateau emerges. Our estimate based on 
data around t = 10 5 suggests that a = 1.08(1). However, 
such an estimation is very close to unity and in addition 
it seems to slowly decrease for increasing t. Such a slow 
decrease suggests that asymptotically a might be equal 
to unity but the power-law scaling is affected by logarith- 
mic corrections. To examine such a possibility, we plotted 
iogioO 2 )/*) vs. bg 10 [log 10 (t)] (Fig.©. The asymptotic 
linearity of the data confirms the presence of the loga- 
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FIG. 4. The time-dependent exponent a(t) calculated in the 
two-dimensional case for several values of TV. 
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FIG. 5. Theplotoflog 10 ((r 2 )/t) vs. log 10 [log 10 (t)] calculated 
in the two-dimensional case. The dashed straight line has a 
unit slope. 



rithmic corrections and, since the slope is close to unity, 
we obtain (r 2 ) ~ tlog 10 (t). 

We used the same technique in the three- and four- 
dimensional versions, and the final results are shown 
in Fig. IS] and Fig. [JJ respectively. In these higher- 
dimensional cases, the linear size L increases slower with 
the number of cities TV, and that is why the curves bent 
downward for smaller values of t. The emerging plateaux 
indicate that a is very close to unity both for d — 3 and 4. 
Let us also notice that in the d = 1 case the salesman, 
possibly after some initial switches, will stick to one side 
and then move ballistically in the same direction, triv- 
ially implying a = 2. 

For the d — 2 and 3 cases, we also analyzed the dis- 
tributions of the length I between successive cities on 
a trajectory. One should take into account, however, 
that these distributions are not entirely stationary and 
to some extent vary in time t (in general, very long dis- 
tances appear in the late-time evolution and they could 
be of the order of the system size L in the final stage) . To 
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FIG. 6. The time-dependent exponent a(t) calculated in the 
three-dimensional case for several values of TV. 




0.5 1 1.5 2 2.5 3 

iog 10 (t) 

FIG. 7. The time-dependent exponent a(t) calculated in the 
four-dimensional case for several values of TV. 
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FIG. 8. The plot of the logarithm' of the distribution P(l) of 
intervals I calculated for paths of the length t = 5000. The 
inset shows that the variability of the average interval (I) (at 
time t) during such runs is rather small. 
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FIG. 9. The same distribution P(l) as in Fig. [8] but in the 
log-log scale. The inset shows the log-log plot of P(l) for runs 
where all cities are visited (N = 10 ). The dotted line has 
the slope -3 (~ l~ 3 ). 



reduce this variability, we generated relatively short tra- 
jectories (most likely none of them reached the boundary 
of the system). 

The logarithmic plot of these distributions (Fig. [8]) 
shows that, in the d = 2 version, the distribution is 
much broader than in the d = 3 version, and perhaps 
of a qualitatively different large-/ tail. Let us notice that 
in both cases (/) is similar and of the order of unity (inset 
in Fig. [5} but the recorded fluctuations of I are several 
times larger in the d = 2 case. Let us also observe that 
large-/ tails of these distributions do not seem to have a 
power-law form. This indicates that the logarithmic cor- 
rections seen in the d = 2 case are not due to sufficiently 
broad (and perhaps Levy-like) distributions of lengths / 
but rather due to the development of long-time correla- 
tions. We will return to this point in the final part of the 
article. 

The log-log plot of these distributions confirms the ab- 
sence of power-law tails for large / (Fig. It turns 
out that the other side of these distributions, namely 
the small-/ tails, can be understood with some simple 
arguments. Indeed, these tails do not differ much from 
the small-/ tails in the initial distribution of cities. For 
simplicity, let us examine first the d = 1 case, where N— 
1 cities are uniformly distributed on the interval of the 
length L — N and one city is placed in the center of 
the interval (i.e., it has the coordinate L/2). Calculat- 
ing the probability density Pi (x) that exactly one city is 
at the distance x from the central city (to the right or 
left) and N-2 cities are further than x from it, in the 
limit N -> oo, we obtain P x {x) = 2(jV L ' 1} (l - ?f)N-2 = 
2exp(— 2x). Generalizing the above calculations to the 
d = 2 and d = 3 cases, we obtain Pi{x) ~ x exp(— 7ra; 2 ) 
and P%{x) ~ x 2 exp(— |7ra; 3 ), respectively. As can be seen 
in Fig. [SI the small-/ tails are indeed close to the tails 
of the initial distributions both for d — 2 and 3. One 



can also notice that for large x the initial distributions 
P%{x) and i-3(ir) decay very fast (faster than exponen- 
tially) while the corresponding distributions P(l) decay 
slower (see Fig. [5]). Thus, the traveling salesman gen- 
erates paths where short distances appear with approxi- 
mately the same distribution as in the initial distribution 
but long distances are generated much more frequently. 

It is perhaps interesting, especially in the context of 
modeling animal search strategies, that the distributions 
P(l) become much broader when calculations are made 
on paths that visit all the cities (inset in Fig. [5]). This 
is because in the final part of the tour the salesman has 
to visit all the "leftovers" that are typically quite distant 
from each other. These numerical data suggest that in 
the two-dimensional data the tail of the resulting distri- 
bution might even have the power-law decay, P(l) ~ Z _Al 
and fj, ~ 3, and thus could be consistent with the limiting 
Levy distribution. Let us emphasize, however, that, for 
paths that visit all the cities, the generated sequence of 
lengths is nonstationary (early- and late-stage statistics 
differ very much), and assembling the data into a sin- 
gle distribution is rather questionable. Let us also note 
that there are some claims that the Levy distributions in 
the foraging data appeared most likely due to a spurious 
inclusion of late-stage flights [12j. These flights, e.g., of 
albatrosses returning to their nests, were obviously very 
long but should not be considered a part of the searching 
pattern. 



IV. FINAL REMARKS 

In this section we briefly summarize the overall behav- 
ior of our model and comment on possible relations with 
some other statistical-mechanics models. The memory 
of the salesman restricts its movements and effectively 
pushes the salesman out of regions that have already 
been visited. For d = 1, this strategy trivially leads 
to the superdiffusive behavior with a = 2. For d = 3 
and 4, the long-time memory turns out to be basically 
irrelevant and an ordinary diffusion (a = 1) is restored. 
The dimension d = 2 seems to be the critical dimension 
separating these two regimes and the logarithmic cor- 
rections are an expected feature in such a case. Let us 
note that a similar memory mechanism generates a su- 
perdiffusive behavior in the so-called self-avoiding walk 
problem [l7| . However, this model of a polymer is known 
to be superdiffusive for d = 1 (a = 2), d = 2 (a = 1.5), 
and d = 3 (a = 1.2), with the normal diffusion resum- 
ing for d > 3. Our model seems to be related to yet 
another model of the self-avoiding walk, introduced by 
Amit et al. [Hj], where the walker tries to avoid lattice 
sites that it has already visited. Renormalization-group 
arguments show that for such a model the critical di- 
mension equals to 2 and the logarithmic corrections are 
present [THEH in the critical dimension. It is likely that, 
at the coarse-grained level, our model is described by the 
same field theory as the model of Amit et al. Its di- 
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rcct applicability to search strategies is rather limited 
but our model might indicate that the very dimension 
of the search space plays an important role. In partic- 
ular, it suggests that much broader distributions should 
appear in the two-dimensional rather than in the three- 
dimensional search spaces. Since a horizontal span of the 
search space for albatrosses or sharks is much larger than 
a vertical one, perhaps this space could be considered as 
effectively two-dimensional and, as our work suggests, 
broad distributions might appear in such a case. Thus 



it would be interesting to examine species for which a 
search space is less anisotropic. 
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